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ABSTRACT 

We present a statistical survey of almost 10 000 radio Type III bursts observed by the Nancay 
Radioheliograph from 1998 to 2008, covering nearly a full solar cycle. In particular, sources sizes, 
positions, and fluxes were examined. We find an east-west asymmetry in source positions which 
could be attributed to a 6±1° eastward tilt of the magnetic field, that source FWHM sizes s 
roughly follow a solar-cycle averaged distribution ^ 14i/ _3 - 3 s -4 arcmin -1 day -1 , and that 
source fluxes closely follow a solar-cycle averaged ~ 0.34 1/ -2 - 9 ^" 1 - 7 sfu -1 day -1 distribution 
(when v is in GHz, s in arcmin, and S u in sfu). Fitting a barometric density profile yields a 
temperature of 0.6 MK, while a solar wind-like (cx h~ 2 ) density profile yields a density of 1.2x 10 6 
cm~ 3 at an altitude of 1 Rs, assuming harmonic emission. Finally, we found that the solar-cycle 
averaged radiated Type III energy could be similar in magnitude to that radiated by nanoflarcs 
via non-thermal bremsstrahlung processes, and we hint at the possibility that escaping electron 
beams might carry as much energy away from the corona as is introduced into it by accelerated 
nanoflare electrons. 

Subject headings: Sun: corona - Sun: radio emission - Sun: particle emission 



1. Introduction 

Particle acceleration events in the quasi-collisionless 
plasma of the solar corona create supra-thermal 
electron beams that propagate along magnetic 
field lines, either away from the Sun into inter- 
planetary space (on so-called "open" field lines), 
or downwards into the chromosphere (towards the 
footpoints of magnetic loops) . Higher-energy elec- 
trons race ahead of the lower energy ones, creating 
a bump-in-tail instability in the particle distribu- 
tion. Landau resonance with the unstable electron 
beams creates Langmuir waves, which are believed 
to undergo non-linear wave-wave interaction and 
produce electromagnetic emissions at the local 
plasma frequency or its harmonic. As the electron 
beam propagates to higher or lower altitudes (at 
spee ds of sa c/3 in the low c orona, to « c/10 at 1 
AU, iPoauerusse et al.l Il99q ) , and hence to lower 



or higher densities and plasma frequencies, there 
is a drift in the frequency of emitted radiation. 
This coherent plasma emission is generally called 
a Type III radio burst. Type Ills have been ob- 
served from frequencies as high as wl GHz at the 
bottom of the corona to 30 kHz at wl AU, and 
even lower further out. Takin g only bursts drifting 
from high to low frequencies, lAlvarez fc Haddock! 
(|l973l ) have fitted bursts in the 74 to 550 MHz 



range by the relation ^ ss -0.01 v l si MHz s" 



Nit a et al 



([20021 ) have investigated the peak 
flux distributions of 40 years of spatially- integrated 
solar radio burst data, from 0.1 GHz to 37 GHz. 
They have found that the peak flux density dis- 
tribution of events, dN/dS„ followed a power-law 
with a negative index ssl.8, similar to that found 
in many X-r ay studies (for rece nt work on the 
topic, see e.g. Hannah et al. 2008h . 



In this work, we have gathered the solar Type 
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(individual ones last about one second). Differ- 
ent groups of Type Ills, even hours apart, are 
however sometimes bundled together in a single 
report. Moreover, certain reported time interval 
are clearly erroneous, probably due to data entry 
error. We have therefore ignored all Type III re- 
ports longer than 10 minutes, (which constituted 
~14% of Type III reports), and were left with 
8931 Type III burst reports from 1998 to 2008. 

The NRH actually produces images at up to 
1/8 s temporal resolution, but these are stored on 
magnetic tapes, and each individual day of data 
must be retrieved manually. The 10-second data 
is stored in a file system, and it is therefore much 
easier and faster to process. This is the reason why 
it is used in this study. This coarse temporal res- 
olution leads to weak and/or short-duration (such 
as second-long Type Ills) bursts being drowned 
in the solar background. Hence, we have decided 
to retain only well-determined Type III bursts for 
this study: those which had peak brightness tem- 
perature of 10 MK or above, well above any quiet 
Sun value. This is the case for 91.4% of cases at 
164 MHz, down to 26.2% of cases at 432 MHz. 
The difference is easily explained by the fact that 
Type III bursts tend to have lower brightness tem- 
peratures at 432 than at 164 MHz (see also ^3.31) . 
A priori, the NGDC-reported burst spectral band 
is not necessarily a good indicator of the actual 
bandwidth of the Type III burst. (NRH, pos- 
sessing «44 dishes, is much more sensitive than 
the single-dish spectrometers that typically report 
to NGDC). On the other hand, we have empir- 
ically observed that taking only bursts either at 
frequencies for which the peak brightness temper- 
ature in the map is above «10 MK, or within the 
NOAA-reported spectral range, lead to very simi- 
lar datasets. 

The 10-second data used in this study also leads 
to cases where several individual bursts are being 
bundled together into a single event. Hence, in 
this work, a "Type III event" is often "a group of 
Type Ills within the same 10-second interval" . 

Frequencies 150.9, 164, 236.6, 327, 410.5, and 
432 MHz (with a 3-dB bandwidth of 0.7 MHz) 
were used by NRH during our period of interest. 
Figure [1] clearly shows a decrease in radio activ- 
ity as the solar cycle approaches minimum, and 

1 ftp://ftp.ngdc.noaa.gov/STP/SOLAR.DATA/SOLAR.RADIO/SPEfepR^B/ in events from 2003 November 5 to 2004 

January 25 is due to NRH being off-line while 



III peak flux densities, source sizes and positions 
at different frequencies, observed by the Nancay 
Radioheliograph (NRH) over the period 1998- Jan- 
01 to 2008-Apr-01, i.e. covering nearly a full solar 
cycle. 

In Section [2] we describe the data selection pro- 
cess and some of its limitations. We describe the 
observations in Section [3] and discuss them in Sec- 
tion 21 Finally, in Section [5j we summarize our 
results and present some conclusions to this study. 

2. Data selection and caveats 

We selected our events from the list of solar 
radio bursts published by the National Oceano- 
graphic and Atmospheric Administration's Na- 
tional Geophysical Data Cente^ (NOAA/NGDC), 
from 1998 January 1 to 2008 April 1, covering 
almost a full solar cycle. We removed all bursts 
that were outside the observing times and spectral 
range (^150 -450 MHz) of the Nancay R adiohelio- 
graph fNRH. iKerdraon fc Delo"uislll997l) . We have 
not gone further than April 2008 because NRH 
started observing the Sun with a different set of 
frequencies, and we wanted a consistent set of fre- 
quencies for our statistical study. The NOAA list 
contains information on individual radio bursts, 
such as reporting station, start and end times and 
frequencies, as well as spectral type and intensity. 
This list is compiled from reports, often gener- 
ated manually, by observers from various ground 
stations around Earth (Culgoora, Ismiran, Lear- 
month, Ondrejov, Palehua, Sagamore Hill, San 
Vito, Bleien, etc.), using spectrometers covering 
varying bands, from as low as ^30 MHz to as 
high as a few GHz. Hence, certain characteris- 
tics such as start and end frequencies, as well as 
burst intensity, are somewhat subject to both in- 
dividual spectrometer sensitivity and individual 
observer's perception. In cases where the same 
event was reported more than once (by different 
observatories), only the first such report was kept, 
and any extraneous ignored. We kept only Type 
III radio burst reports within the times («8:30 to 
15:30 UT, varying during the year) and the spec- 
tral band (150-432 MHz) where the NRH was 
observing the Sun. Groups of decimetric Type 
Ills do not last much longer than a few minutes 
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Fig. 1. — Black: Time series of the quarterly number of NOAA-reported solar radio bursts, at times and 
frequencies when and where NRH was observing. Red: Sun Spot Number. Blue: F10.7 index (in sfu). 
Three-month time bins or smoothing window were used for all displayed data. NRH was not observing at 
150.9 MHz before mid-2002, and was not operating from November 2003 to January 2004 (see text in $2]for 
more details). 



anti-alias antennas were being added to the ar- 
ray. During our decade of interest, the 150.9 MHz 
frequency has gradually replaced the nearby 164 
MHz, no longer reserved under French law. Hence 
the observations at 150.9 MHz have started later, 
with an overlap of a few years with the (now no 
longer used) 164 MHz. It can be noted that the 
oc currence rates mat ch qualitatively that 
of lLobzin et al.l (|201ll ). including an almost 
^two-year periodicity, unexplained so far. 
Using NRH 10-second images, and a fixed 30" 
pixel resolution at available frequencies (this stan- 
dard coarse resolution, still well below the instru- 
mental resolution at 432 MHz, was chosen for ease 
of data handling), the following were determined 
(among other things): time of peak flux density, 



location and brightness temperature of brightest 
pixel in the source, and 2D gaussian fitting pa- 
rameters to sources present in the image. From 
the latter, source fluxes can be deduced by com- 
puting the gaussian volume. It is important to 
note that each frequency was treated mostly inde- 
pendently from each other, see §2.11 and 12.21 for a 
discussion. 

Paragraphs 2.1 to 2.6 detail some caveats in our 
selection and data analysis process. 

2.1. Time interval of reports: 

Firstly, as stated earlier, there can be several 
Type Ills temporarily close to one another, within 
the 10-second accumulation of the data at hand, 
leading to confusion and averaging of their charac- 
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teristics. Secondly, our methods looks only for the 
strongest burst within the NOAA report's time in- 
terval (there can be several inside the time range 
reported), and does so for each frequency inde- 
pendently. This has for consequence that we are 
missing any other bursts within the reported time 
interval, and that the burst characteristics derived 
at different frequencies might not relate to the ex- 
act same Type III burst. About ^50% of radio 
bursts from the same report have hence a time 
difference greater than 10 s between 164 and 432 
MHz). This has an obvious effect on the abso- 
lute occurrence rates of Type III bursts (they are 
probably higher than reported here), but we esti- 
mate the influence on the power-law shape of their 
distributions (see $3]) to be minimal. 

2.2. Alias ambiguity: 

For bursts occurring before the installation of 
anti-alias antennas at NRH (between 2003 Novem- 
ber 5 and 2004 January 25), there can be an ambi- 
guity in source positions at frequencies above 164 
MHz. This is because an alias is sometime present 
in the reconstructed NRH image, and sometimes 
has the brightest pixel in the map (particularly 
when the real source happens to be beyond the 
solar limb, but the alias lies on top of the solar 
disk). The distance between the real source and 
an alias decreases with increasing observing fre- 
quency, and any alias is beyond the imaged field 
of view (and hence no alias ambiguity exists) at 
frequencies <164 MHz. It is usually easy to dis- 
tinguish between a solar source lasting at least a 
few minutes and its alias: the alias moves quickly 
(with time) in a straight line across the solar map. 
But this method is of little use in our study, as 
we examine events that seldom last longer than a 
single 10-second frame. We will instead use the 
fact that the true source is likely to be spatially 
near the sources at other frequencies, assuming the 
radio bursts spans the necessary frequency range. 
For maps taken at progressively higher frequen- 
cies than 164 MHz, we determine whether there is 
possiblity of an alias in the map in the following 
manner: after locating the source with the bright- 
est pixel, we determine whether there is another 
source in the map at an appropriate alias posi- 
tion, and with similar peak intensity (within 20%, 
which corresponds to a 2 MK difference for the 
minimal 10 MK burst). The "correct" source is 



deemed to be the one whose position is closest to 
the position of the "correct" source in the map 
of the next lower frequency. This simple method 
is not perfect: in some cases, emission at lower 
frequencies is simply not there, because the con- 
ditions for plasma emission and/or propagation 
to Earth are not adequate, and hence, at higher 
frequencies, the method may home in on other 
sources. 

Before the installation of anti-alias antenna, 
there was a potential alias confusion (two sources 
of similar intensities, at correct alias positions, as 
described above) in our data 58% of the time at 
432 MHz, progressively lower at lower frequencies, 
down to 0% at 164 MHz and below. 

In summary, source locations before January 
2004, low frequencies (150.9 and 164 MHz) should 
be always accurate, while higher frequencies may 
have errors in some of the positions, but we expect 
the impact on the statistical distributions of burst 
locations ( ^3.21) to be minimal (see coronal den- 
sity fittings in $4] for a confirmation). The study 
of source sizes ( £!3.1j) and fluxes ( EI3.3j) should be 
unaffected by this issue, as is all the data taken 
after January 2004. 

2.3. Bursts fluxes: 

The so-called "Quiet Sun" flux is of the order 
of 10 sfu at NRH frequencies (and brightness tem- 
peratures at Sun center of ~0.7 MK at 164 MHz). 
Hence, a lone, 1-sfu, 1-second long Type III burst 
with area about 1/100 the solar disc (i.e. smeared 
out to about the size of the NRH beam at 164 
MHz, about 3.2' FWHM) would produce an 
additional surface brightness equal to that of the 
quiet Sun in a 10-second image. Successfully fit- 
ting a gaussian to such a weak burst becomes dif- 
ficult, and explains the turnover at low fluxes in 
the distributions presented in i)3.3l 

2.4. Astrophysical sources: 

The Crab Nebula, or one of its alias, is the 
most obvious astrophysical source which can enter 
the NRH field of view (around June each year). 
Its surface brighteness is similar to the Sun at 
NRH frequencies, i.e. much lower than the 10 MK 
threshold used for our selection. Hence, we expect 
such sources to have little impact in our study. 
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2.5. Refractive effects: 

Earth ionospheric refraction can significantly 
perturb the apparent position of a source in 
ground-based observatories, sometimes exceeding 
several arcmin s in the 100-3 00 MHz range (e.g. 
lBougeredll98lt lMercier|[l986h . These effects are 
mostly random position shifts due to ionospheric 
gravity waves. They are not correlated with the 
positions of sources on the Sun, and should not 
affect our statistics. We have therefore ignored 
those effects in the present study. 

3. Data Analysis 

We have studied in detail the distribution in 
spatial sizes f£ )3.1[) . in positions f ^3.2p . in fluxes 
(33T3J) , and the rank-order cross-correlations be- 
tween several of the observables ( 33. 4p . 

3.1. Burst spatial size distribution: 

Radio sources were fitted with 2D elliptical 
gaussians. These yield observed gaussian sizes 
o~ a ,obs and Ub.obs (semi-major and semi-minor e- 
folding lengths), and a tilt angle o b s of the semi- 
major axis with respect to the map x-axis (so- 
lar east-west). The FWHM values, s a .obs and 
Sb,obs, are obtained by multiplication of <7 a}0 bs and 
o~b,obs with 2\/2 In 2 ps 2.355. In Figure [2j we 
plot a histogram of the rms averages s rms . b s = 

\J s a obs + S l obs- The roll-over at low source sizes 
is due to the interferometer beamwidth, which is 
about 3.2-5.5' FWHM at 164 MHz, proportionaly 
smaller at higher frequencies, and varying with the 
season and time of day. 

Assuming that the observed source, the inter- 
ferometer beam, and the true (deconvolved) source 
were all elliptical gaussians (with different tilt an- 
gles), we have recovered the true source size (see 
Appendix |A1 for mathematical details), and plot- 
ted the results in Figure [3] and their averages in 
Table [TJ Note that some of the largest events (in 
size) are probably multiple sources, but it was not 
the purpose of this statistical work to discriminate 
between such. 

As the spectral indices of the distributions for 
the bursts selected in Section [5] were very close 
(from 4.4 to 4.9) at all frequencies, we decided 
to plot the normalization constants A v as a func- 
tion of frequency. But to remove any normaliza- 
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Fig. 2. — Histogram of the radio sources' rms 
FWHM s — s rmS}0 bs, at different frequencies. The 
red curves are fittings =p = As a to the the his- 
tograms, from the bin to the immediate right of 
the peak in the histogram and above. The blue bar 
delimits the minimum NRH point spread function 
(which can vary by almost a factor 2, depending 
on the season and time of day). 
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Fig. 3. — As in Figure [2] but displaying true 
(deconvolved) FWHM rms source sizes. The red 
line and labels are power-law fittings to the data. 
The blue bar delimits the minimum NRH point 
spread function (which can vary by almost a factor 
2, depending on the season and time of day). 



tion issue stemming from the fact that not all fre- 
quencies were uniformly employed by NRH during 
the solar cyle (particularly, 150.9 MHz was used 
only after 2002/10/28), we have used a slightly 
smaller subset of our data to which we have fit- 
ted similar power-laws: we have taken only those 
events that occurred when NRH was observing at 
all six frequencies simultaneously, and plotted in 
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Table 1: Mean and standard deviation of observed 
and deconvolved (or "true") source sizes: 
Frequency: Observed Deconvolved 
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size ['] 


size ['] 
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Fig. 4. — Scatter plot of the constants A obtained 
from fittings similar to that of Figure [3J 



Figure |4] the normalization constants A v of these 
new power-law fits. Figured] suggests a power-law 
frequency dependence of the normalization con- 
stant Am B , which has led us to infer that the 
distribution for true source FWHM size s could be 
of the form: 

<'V 



ds 



= Bv' L 



We have done two-dimensional fittings of our data, 
using Equation [T] as model. The fitting parame- 
ters (B, a, (3) have been found to be moderately 
dependent on the choice of fitting intervals (and, 
to a lesser extend, the binning). We therefore 
prefer to give a range of values for which \ 2 
near unity was obtained: (B, a, 0) from (2.0,-3.0,- 
2.8) to (20.0,-5.0,-3.8), for an average of «(7,-4,- 
3.3) over the 2002-2008 period. The result is in 
arcmin" 1 day -1 if v is in GHz and s in arcmin. 
To account for the factor ss2 ratio between the 
1998-2008 and the 2002-2008 average burst rates 
(Figure [T]) , the normalization constant should be 
changed from 7 to 14 to get a more proper solar- 



cycle average. 

When the Sun is low on the horizon during cer- 
tain period of the year, the beam shape can be ex- 
tremely asymmetric. To insure that these extreme 
cases have not influenced our results, we have 
therefore run the same study, removing all events 
occurring in November, December, and January, 
as well as all events within two hours of sunset or 
sunrise during February and October. The differ- 
ences with the above results were negligible. 

A more thorough study of source shape and 
structure, which would address e.g. c ases where 
more than a single source is present (e.g. lPick et al. 
19981 . and references therein), is beyond the scope 
of this paper. 

3.2. Radioburst location: 
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Fig. 5. — Occurrence density distributions for the 
radio bursts selected in Section H on a 1' pixel 
grid. The number in upper right corner is the to- 
tal number of events with reliable positions. The 
intensity scale is linear from white to black, with 
black representing the maximum pixel value as 
specified in each plot. The arrows indicate the 
±15 degrees heliographic latitudes. 

Figure [5] shows the distribution of radio bursts 
for our selected bursts, using the positions that are 
likeliest to be correct (as discussed in Section [2]), 
while Figures El and |8] show the same for three 
different characteristic years in the solar cycle. 

The most striking feature is, at high frequen- 
cies, the concentration of radio bursts in two dis- 
tinct bands (around latitudes ±15 to ±30 de- 
grees for 432 MHz) during early solar maximum 
years (Figure [B]) , and at lower latitudes as the so- 
lar cycle advances. This is in remarkable agree- 
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Fig. 6. — As Fig. taking only events that oc- 
curred in 1998. No observations at 150.9 MHz 
were made that year. 
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Fig. 7. — As Fig. [5j taking only events that oc- 
curred in 2002. 
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Fig. 8. — As Fig. taking only events that oc- 
curred in 2006. 



ment with observed active region and microflarc 
positions and their solar-cycle dependence (e.g. 



Christe et al.l l2008al ). Later in the solar cycle, 



bursts are located mostly south of the solar equa- 



tor, just as active regions were (Higgins et al., 
2012, in preparation). 

The second most striking feature is the system- 
atic shift of about 2' =120" westward of the mean 
position of all radio bursts (Figure|9]). While a sys- 
tematic shift of instrumental origin of up to 0.3' 
westward for all frequencies was expected (from 
comparison between VLA and NRH positions of a 
radio spike - priv. comm. P. Grigis.), clearly it 
cannot account for all of the observed shift. 
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Fig. 9. — Spatial distribution of all radio burst 
positions between 1998 and 2008: histograms for 
all six frequencies (top six: east- west coordinate, 
bottom six: north-south coordinate) with means 
and their errors. A westward shift between 115 
and 149" can be observed at different frequencies. 



3.3. Source peak brightness temperatures 
and fluxes 

In this section, we investigate the distribution 
of "peak" (in terms of 10-second averages) bright- 
ness temperature Tg (in K) and source fluxes S u 
(in sfu, where 1 sfu = 10~ 22 W m~ 2 Hz -1 ). Fig- 
ures [10] and QT] show that the peak brightness tem- 
peratures and fluxes of our radio bursts, follow a 
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law, very close to what iNita et al 



found Equivalently, the number of bursts 
above Tb or S v follow a power-law with negative 
index one higher (see Figure [T2l for S v ). The roll- 
over at low fluxes is most likely due to our selection 
criteria (see discussion in Sj2]). 
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Fig. 10. — Histogram of deconvolved (i.e. assum- 
ing both observed and true sources are gaussian- 
shaped) peak brightness temperature of all Type 
III bursts. The dashed blue line is a power-law 
■ ^L = AT% fit to the data, using the C-statistic 



Cashlll97<i . technically better suited than Pois- 
son statistics for datasets with small number of 
counts per bin (as is our case for the high-value 
bins, but in this case leading to negligible differ- 
ences). The associated best-fit parameters are in 
the lower left corner. 



Note that unresolved sources always have lower 
brightness temperature than their "real bright- 
ness" temperature, and that a source's "real" 
brightness temperature is always smaller or equal 
to its true temperature (opacity). For these rea- 
sons, the flux S u is probably a less misleading 
quantity to use than Tb, and we will concentrate 
on S v in the remainder of this section. 

As the spectral indices in Figure [TJJ at differ- 
ent frequencies were very close, we decided to plot 
the normalization constants A v as a function of 
frequency (Figure [T5|) . As with the source sizes in 
H3.ll we have taken a subset of our data, using only 
events occurring when NRH was observing with 
all six frequencies. Power-law fittings to these new 
distributions yield more meaningful normalization 
constants A„, plotted in Figure [TJJJ Again, notice 
that the coefficients A v appear to have a power- 
law dependence with frequency: A v w Bv$ . We 
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Fig. 11. — Histogram of gaussian source fluxes 
(averaged over 10 seconds), with power-la w fit- 
tings 4jp = AS" using the C-statistic ( Cash 
1979 i. 
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Fig. 12.— Cumulated data N(> S u ). The blue 
dotted line is the integral of the fitting found in 
Figure QTJ while the green dashed line and coef- 
ficients ar e derived from the Ma ximum-likelihood 
method o f C rawford et al.l (|l970l ) . They match al- 
most exactly. 



have therefore attempted to fit all our data with 
the following model: 



dN 



B i/* 3 



(2) 



The best-fitting parameters for the 2002-2008 pe- 
riod were found to be (with v expressed in GHz, 
and S v in sfu): B=0.17±0.01 sfir 1 day" 1 , a=- 
1.7±0.05, /3=-2.9±0.1. These fitting parameters 
have proven to be very stable (and all with near 
unity \ 2 ) even if changing fitting intervals, and 
hence much more robust than those found in Sec- 
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Fig. 13. — Plot of the normalization coefficients 
A for fittings similar to Figure |T2] but taken only 
over times when all six frequencies were observing, 
as a function of frequency. 

tion 13.11 for size distributions. To account for the 
factor «2 ratio between the 1998-2008 and the 
2002-2008 average burst rates (Figure Q]), the nor- 
malization constant should be changed from 0.17 
to 0.34 to get a proper solar-cycle average. 

It is well known that more Type III bursts are 
seen at low frequencies than at high frequencies, 
and the above relat ionship reflects this. In fact, 
Dul k et al.l (|200lh found the spectrum of a 
Type III burst in the 3-50 MHz range to 
have a negative spectral index close to 3, in 
good agreement with our results. 

3.4. Rank-order cross-correlations: 

We cross-correlated most of the observables 
mined from our study. The results and method 
of correlation are described in Appendix |EJ The 
only noteworthy relationship that we have ob- 
served is the anti-correlation between the "radial 
elongation" (the ratio of the deconvolved source 
size in the radial direction over the deconvolved 
rms source size) of sources and the radial offset 
from Sun center (Figure [T4]) . This anti-correlation 
means that sources tend to be circular near disc 
center, and elliptical near the limb, with minor 
axis along the radial direction. We will discuss 
the implications of this mild correlation in the next 
Section. 

There is a noteworthy absence of strong corre- 
lation between the source flux S v and the source 
size strue, which will also be discussed in Section[4] 




0.0 0.5 1.0 1.5 2.0 

Eost-Wesl distance from Sun center [Rs] 

Fig. 14. — Scatter plot of radial distance vs. "ra- 
dial elongation" (ratio of deconvolved source size 
in the radial direction over deconvolved rms source 
size) of sources. The mild anti-correlation indi- 
cates that, the further away from Sun center, the 
more elongated in the direction perpendicular to 
the radial sources tend to be. 



Due to the longer path in the corona that emis- 
sions from sources near the limb have than those 
near disk center, it is expected that turbulent scat- 
tering would make them appear slightly larger 
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on av erage than their on-disk brethren (jBastian 
1994). There appears to be indeed a slight cor- 
relation («0.2) between deconvolved source size 
Strue and angular distance r from Sun center in 
our data (particularly at the lower frequencies). 

4. Discussion 

Source sizes: The variation of source sizes with 
frequency relies on a convolution of different ef- 
fects: the magnetic field structure opening as a 
function of height in the corona, the fact that we 
may have several different Type Ills at slightly 
different locations in the same 10-second time in- 
terval, the distribution of intrinsic sp atial width 
of the electron beam (jBastian 1994 . considered 
point sources in his work), and, of course, a host 
of (most ly refractive in origin) propagation effects 
(see e.g. lPoquerusse fe McIntosblll995l:[Dulkll2000l 
for recent overviews). For example, true source 
sizes could be 3.5 (harmonic) to 1.7 (fu ndamen- 
tal) s maller than observed, according to iMelrosei 
(Il989h . 

Average values of Type III burst sizes have been 
measured at frequencies below 169 MHz in the 
1970s-1980s dBqugeret et al lll970t [steward 1 1974; 



Dulk et alJ[l979HDulk fe Suzukilll980h . Thev are 
found to increase with decreasing frequencies and, 
on the average, to be of the order of 5 ' at 169 MHz 



(valu e s between 2' and 7 \ see also IZlobec et al 



1992 : Mercier et al. 20061 ). which compares well 
with the 4.5±1.6' value in Table [TJ The first 
spatially resolved observations of Type III bursts 
at high time resolution with the Nangay Radio- 
heliograph showed that Type III bursts can be 
resolv ed into narrow components ( Raoult fe Pickl 
1980j) (elementary size of the order of 2') with 
a typical size which increases with time during 
the burst, potentially reaching 5' to T at then 
end of the burst. What we have measured in the 
present paper where we do not use high time res- 
olution data will reflect more the total size of the 
whole Type III burst group rather than the size 
of individual components. Systematic studies of 
burst sizes have been performed only at 164 MHz 
(|Raoult fc Pickl Il980t IPick fc Jilll987h . The dis- 
tribution of east-west sizes showed that Type III 
bursts may be resolved in a few cases in compo- 
nents with real sizes around 2', but that the dis- 
tribution of sizes can reach 9'. The present study 



provides the first information on the sizes of Type 
III bursts at frequencies above 169 MHz. It shows 
that Type III bursts are smaller at high frequen- 
cies than at low frequencies and that the size de- 
creases as v~ 
Theory (e.g 



-3.3 



Bastianl 11994 and references 



therein) predicts that scattering of the emission in 
the solar atmosphere would lead to a v~ 2 depen- 
dence of the source size. This has been confirmed 
by observa tions of cosmic sources through the so - 
lar corona (|EricksorJl96ilColes fc Harmonlll989l) . 
However, previous observations have showed a 
quasi-linear dependency of Type III source sizes 
(from saO.l to »150 MHz) wi th wavelength 
(jSteinberg et all Il985l iDulkl l2000h . Our obser- 
vations support neither. 

Although he himself states his analysis becomes 
inva lid for radio w avelengths longer than a few 
dm, Bastianl ( 1994 ) predicts a sizeable increase of 
source sizes near the limb. The fact that we seem 
to observe only little of it might be an indication 
that there is a cutoff in the level of turbulence. I.e. 
that below a certain altitude, turbulence would no 
longer be proportional to the ambient density as 
it is for distances greater than 1.7 Rs- 

Weak correlation between source sizes and 
fluxes: This is consistent with the interpre- 
tation that emission comes from a spatially 
smal l region and that propagatio n effects (see 



e.g. iPoquerusse fc Mcintosh! Il995l for a short 



overview) are responsible for the lack of good cor- 
relation. Notice, however, that the correlation 
coefficient gradually improves at higher frequen- 
cies (from 0.08 at 164 MHz to 0.35 at 432 MHz), 
suggesting that these (refractive) propagation ef- 
fects become less important at higher frequencies. 

Source anisotropy near the limb: The slight 
anisotropy observed as sources get nearer to the 
limb can be interpreted as being due to the fact 
that they are being observed more edge-on than 
when at Sun center, and that the emitting layer is 
thinner than the horizontal source size. 

Source location: Type HI emission at the lo- 
cal plasma frequency (fundamental) is expected 
to be radiated primarily in a dipolar pattern, 
perpendicular to the local magnetic field. For 
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radiation at the harmonic, a quadrupolar pat- 
tern is expected. However, the propagation 
of radio waves is controlled both by large and 
small scale structures which can modify the pri- 
mary directivity of the radio emission thr ough re- 
spect ively refraction (including ducting, iDuncan 
1979) and scattering processes. Particularly, 
emission at the fundamental tends to become 
aligned with the local d e nsity g radient. S ee e.g . 
Zheleznyakov fc Zaitsevl (Il970l): [Melrose! (|l986l ); 
Cairns! dl987allbl): iRobinson et al.l dl994l) for theo- 
retical work, IPoq uerusse fc Mcintosh! (119951) for a 
recent short overview, and lTheiappa et all ( 2012 ) 
for recent observational work. 

The first investigation of the directivity of Type 
III bu r sts was achieved by Caroubalos fc Steinberg 
(1974); Caroubalos et al. ( 1974) using stereoscopic 
measurements at 169 MHz obtained simultane- 
ously with the STEREO- 1 experiment and with 
the NRH. Some directivity was found for the Type 
III bursts (especially for the fundamental emis- 
sion). The directivity was however found to be 
less than for Type I bursts. Stereoscopic measure- 
ments of the directivity have then been performed 
using ground-based measurements at 150 MHZ 
and measur ements on ULYSSES in the 1.25-940 
kHz range (jPoauerusse et al.l Il996t iHoang et al 



1997). They found that the average pattern of the 
Type III bursts at low frequencies is shifted 40° 
eastward of the radial direction. The shift and 
width is found to d ecrease with low frequencies. 
Bonnin et al.l ([20081 ) further investigated the di- 
rectivity of Type III bursts in the interplanetary 
medium using calibrated Wind-Ulysses observa- 
tions in the same frequency range. They con- 
firmed an eastward shift of 23° at 940-740 kHz 
and 55° at 55-104 kHz, i.e. increasing with wave- 
length. In all these papers, the shift is attributed 
to a transverse density gradient created by the 
fast wind (propagating along spiraled open field 
lines) overtaking the (mostly radial) slow wind. 
At meter wavelengths, an eas t-west a s ymmetry fo r 



noise storms w a s reported b y i Fokkerl (il960lll963lh 
Suzuki (1961); Le Squerenl ( 1963 ). (see ElgarOv 
19771 pp.111-115. and references within for a re- 
view). The asymmetry is about 0.1 Rs towards 
the west. One of the possible explanations which 
was proposed is that there is a tilt towards the 
east in the magnetic structures. 

In the present study performed at decimet- 



ric/metric wavelengths, we do not observe a oc A 
behaviour of the westward shift in our dataset, but 
this could be due to the fact that NRH observes 
bursts at much higher frequencies and much closer 
to the Sun. On the other hand, combining the ob- 
served ~2' westward shift with a simple geomet- 
ric model (details in Appendices [C]), an average 
tilt angle between the direction of emission and 
the radial to Sun center can be derived, and has 
been found to be around 6 degrees. (This result 
is consistent with a more elaborate model using a 
Monte-Carlo approach, detailed in Appendix ID1) 
This is in the same order of magnitude as was 
found from noise storms asymmetries. If one as- 
sumes that the optimal direction of emission sta- 
tistically corresponds to the direction of the local 
magnetic field, one hence concludes that the mag- 
netic field is on average tilted eastward by «6° 
with respect to the radial direction, at altitudes of 
a few tenths of a solar radius, where emission at 
150-432 MHz occurs. 

Derivation of coronal density profile: In the 

following, we assume that emission is either at the 
fundamental or the harmonic of the local plasma 
frequency, though the latter is generally thought 
to be the dominant emission mechanism in the 
case of decimetric Type III emission. Using posi- 
tion data after 2004/01/25 (the best position data, 
obtained after the installation of anti-alias anten- 
nas), it is possible to derive statistically the aver- 
age height difference between e mission at different 



freque ncies. See Appendix A of lSaint-Hilaire et al 



( 2010t ) for comprehensive details of the method. 
This method assumes that sources at different fre- 
quencies tend to be radially distributed at different 
heights. The effect of the small «6° tilt angle is 
neglected. In the first five panels of Figure [TSJ we 
have plotted all average R (plane-of-sky, or POS, 
distance from Sun center) vs. AR (POS distance 
between centroid at two different frequencies) for 
emission at 150.9 MHz versus the other five fre- 
quencies. The fitted slopes are proportional to 
the average height difference (the nearer to disk 
center, the smaller the height difference, and the 
nearer to the limb; the larger the height differ- 
ence). Note that the line must go through the 
origin, making the fittings much stronger than ini- 
tially appears. The height differences are plotted 
in the last panel of Figure [151 We have fitted two 
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Fig. 15. — Using the methodology presented i n 
detail in Appendix A of lSaint-Hilaire et al. ( 2010l ). 
we have statistically derived the relative heights 
between emission at different frequencies (first five 
plots). In the bottom plot, we have plotted the 
height differences, with error bars, and fitted sim- 
ple density profiles. See text for details. 



enced by whether the emission is at the fun- 
damental or harmonic of the local plasma 
frequency. Best fitting yields 77=31.5 Mm 
(the red curve in the bottom plot of Fig- 
ure [T5]). This scale height corresponds to a 
«0.6 MK corona, lower than the more usual 
quiet Sun values of 1-1.5 MK. This could 
be due to the fact that we are observing 
Type Ills near regions of open field lines, 
and/or the fact that we are near altitudes 
wher e solar wind-like con ditions start to take 
over. iDavid et all (|l998t ) indeed found very 
low temperatures at similar heights in polar 
coronal holes. 

2. A "solar wind-lik e" atmosphere n (h) 



C{h/R s )-' 2 model (jCairns et aU2009l ). which, 
in order to be used on our data, is trans- 
formed into: 



Ah i:j = hi-hj = 8980 mRsVC ( — - — 



(4) 

Where m=l or 2 for emission either at fun- 
damental of harmonic. The normalization 
constant C can be viewed as the electron 
density at an altitude of 1 R$- Best fitting 
yields C=5xl0 6 cm -3 (emission at funda- 
mental) and C=1.2xl0 6 cm -3 (emission at 
harmonic), and the blue curve in the bot- 
tom plot of Figure [15J The numbers for 
harmonic emission are comparatively low, 
while for fundamental emission, they are 
comparable to the Baumbach-Alle n formula 
and o t her published work : e.g 



Ncwkirk 



(|l96ll ); iMann et al.l (Il999h for the o retica l 



mod els, and see Fainberg fc Stone] (1974) 
and Aschwanden fc Actonl ( 20011 ) for a re- 
cent compilation of observational work. 



simple density profiles to this data: 

1. A hydrostatic exponential atmospheric n = 
nae~ h / H model (where H is the scale height, 
and h the altitude above the photosphere), 
which, in order to be used on our data, is 
transformed into: 



Note that we have also used the pre-2004 data, 
and have found there was very little difference with 
the "better" post-2004 data. Both fittings appear 
equally good (or equally bad), and this method- 
ology cannot decisively choose which of these two 
models is the most appropriate. 



Ah i:j = /i 4 



hj = -2H]n— (3) 



Where i,j refer to different frequencies of ob- 
servations. Note that this model is not influ- 



Source intensities and fluxes: Maximum ob- 
served brightness temperature Tb in interplane- 
tary Type III bursts are reported to be sslO 15 K, 
and up to «10 12 K for coronal Type III bursts 
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(lMelrosdll989h . We observe a few events >10 11 K 
in our dataset. 

We have found a power-law dependence of the 
occurrence rate of Type III bursts of peak flux S v 
with a negative s p ectral index of 1.7, very sim- 



ilar to iNita et al.l ( 2002 ) , and to what is found 



in X-rays for flare e nergies (ICrosbv et al 



Hannah et ai1l2008l ). lEastwood et al.l (|2010l ) have 



1993; 



found a 2.1 index for interplanetary Type III so- 
lar radio storm in the 0.125-16 MHz range.. A 
1.7 negative power-law index is also an oft-seen 
value in self-organized c riticalitity (SOC ) stud - 



ics 



■Baket all (Il988l): ILu fe Hamilton! 
Vlaho s et al.l ~(ll995fl : iGeorgoulis fc Vlahosl 



1991) 



1996) 



Aschwandenl (120121 ) . In their study of Type I 



(2012 
, Merci 



noise storms, iMercier fc Trotted (|1997l ) had found 
negative power-law spectral indices of about 3 
to 3.5 for their and have attributed it to 

the predicted low- energy part of SOC theory 
(jVlahos et al.lll995l ). (Type I noise storm bursts 
are indeed less energetic but much more frequent 
than the decimetric Type III bu rsts studied here.) 
Similarly, iMorioka etabl (|2007l ) found a -3.6 in- 
dex for interplanetary (wMHz) "micro-type-III" 
bursts. 

Additionally, we find a cx v& (with j3 «2.9) 
spectral dependence of -f^- , reflecting the fact that 
Type Ills are more easily observable, or occur 
more often, at 164 MHz than at 432 MHz. This ex- 
ponent could provide c onstraints to Type I I I emis - 
sion models su ch as in lRobinson fc Cairnsl ([19981 ): 
Melrose! (|l989l) . 



Solar-cycle averaged radiative output: To 

put things in perspective, we compared the global 
energy radiated by the Type III bursts and the en- 
ergy input of flares into the corona. We have plot- 
ted in Figure [16] the "Radio Burst Occurrence 
Rate" (RBOR) of our observed Type Ills and 
compared it to the energy input from EUV/X-ray 
eve nts to the coro na, derived in other studies (see 
e.g. lHudsonlll99ll for a discussion on this topic). 
The RBOR can be estimated from the power- 
law fitting parameters by using the at 164 
MHz (Figure ITTj). by assuming a burst band width 
of « 300 MHz (|lsliker fc Bend Il994l ). tak- 
ing At=10 seconds (since all our numbers are 10- 
second averages), and computing the total energy 
radiated as observed from D=l AU. We have fur- 
ther normalized to the photospheric surface area, 




10 10 z0 10 zz 10 z 10 z5 10 ZB 10 30 10 3 ' 



Energy [erg] 

Fig. 16. — Black/gray: Approximate flare heating 
input to the solar corona as derived from various 
EUV/X-r ay nanoflare and mic roflare studies: RH: 
RHESSI [ Hannah et aD (<2008h ; SS: Yohkoh SXT, 
IShimizul (119951); TA: TR ACE, IXschwanden et al 
( 20001) ; TP: TR ACE. iParnell fc Jupp| (l2000h 



Hannah et al 



2008) 



EB: SOHO EIT iBenz fc Kruckerl (|2002l) (after 



Red lines: Radiated ener- 
gies. The dotted red line corresponds to the ra- 
diated energy in non-thermal HXR that roughly 
accompanied the RHESSI-derived flare input, us- 
ing an approximated nonthermal bremsstrahlung 
efficiency of 10~ 5 . The solid red line is the radi- 
ated energy derived from this Type III study. 



in order to compare to the "Flare Frequency" of- 
ten used in coronal heating studies: 



RBOR 



1 



1 



dN 



A-kRI AitD 2 AvAt dS v 



(5) 



(Note that A-kD 2 AvAtdS^ = dE is the radiated 
energy for one event). The RBOR is hence the 
occurrence rate of radiated energy distribution 
of Type III events, per unit (photospheric) so- 
lar area. One notices that the radiated energy 
by decimetric Type Ills is about five to six or- 
ders of magnitude less than the energy input by 
EUV nanoflares in the corona (we loosely refer 
to flares observed in EUV only as "nanoflares" , 
and those observed also in X-rays as "microflares" 
or "flares"), about the same ratio that exists be- 
tween the non-thermal hard X-ray (HXR) emis- 
sion associated with X-ray flares and microflares, 
though the emission mechanism is believed to 
be vastly differen t (coherent plasma e mission vs. 
bremsstrahlung). iRamesh et al.l (|20101 ) have esti- 
mated that a 40 sfu metric (77 MHz) interplane- 
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tary Type III burst was produced by a « 3 x 10 24 
ergs electron beam. Using the same assump- 
tions on Av, At, and D as before, this would 
imply an efficiency of non-thermal radiated en- 
ergy over beam energy of «10~ 5 (within an or- 
der of magnitude). It is interesting to note that 
this ratio is somewhat similar to the one be- 
tween radiated non -thermal X-r ays and the en- 
ergy in thick-target ( Brownlll971 ) accelerated elec- 
trons (see Appendix |B]|, To our present knowl- 
edge, there is no obvious reason for the non- 
thermal bremsstrahlung efficiency (an incoherent 
process) and the efficiency of Type III plasma 
emission (a collective process) to be similar, and 
it may be a coincidence. Clearly, more studies 
are needed to acertain radiative efficiencies, but 
assuming they are indeed similar, this has for con- 
sequence that nanoflare and possibly microflarc 
electrons could be similar in number and ener- 
gies as Type Ill-producing escaping coronal elec- 
tron beams! This is in contradiction with recent 



studies ( Christe et al.l 2008b; Saint-Hil aire et al 



2009), which have shown that there is probably 
a factor ss500 more flare electrons than escaping 
electrons. However, note that these studies were 
centered on one-to-one association of precipitat- 
ing flare electron beams and escaping interplane- 
tary electron beams (loosely assuming that inter- 
planetary beams and their counterparts in the low 
corona have similar properties), while the present 
study is more global (and, for example, it allows 
for temporal displacement between the two phe- 
nomena and different rates of occurrence). One of 
the implications is that escaping electron beams 
carry energy away from the corona, partially or 
wholly negating coronal heating by nanoflares. We 
dare not speculate any further without a bet- 
ter handle on radiative efficiencies of Type III- 
producin g electron beams. Th e upcoming NuStar 
mission ( Harrison et al. 20051 ). and perhaps also 
the FOXSI rocket flight, should reveal how many 
electrons are in Type III beams (via their X-ray 
emission), yielding the efficiency of Type III radio 
emission. 

5. Summary, conclusions, and future work 

We have studied solar Type III radio bursts ob- 
served by NRH during 1998-2008, almost a full so- 
lar cycle. From this very rich dataset, we have de- 
termined source sizes, locations, fluxes, and their 



frequency distributions. 

The source size distribution ^can be fitted 
with a double power-law in v and in s. The 
Type III burst sizes are found to decrease as i/~ 3 3 , 
which could reflect a combination of the magnetic 
field opening as a function of height in the corona, 
and the distribution of the spatial widths of elec- 
tron beams in the corona. 

The source flux distribution follows a -1.7 
power-law at all frequencies. A two-dimensional 
power-law fitting yields oc is~ 2 9 S^ 1 ' 7 . These 
values offer additional constraints to theoretical 
models of Type III emission. 

Type III emission generally come from active 
regions (just as flares), and we have observed 
an east-west asymmetry in their location, which 
could be explained by an eastward tilt of the 
magnetic field (compared to the local radial to 
suncenter) at the typical height where our met- 
ric/decimetric Type III emission occurred. We 
have estimated that the tilt angle is about «6° 
at altitudes «0.3-0.6 i?s- 

A barometric fit to the data led to a ^32 Mm 
scale height, and temperatures similar to what can 
be found in coronal holes, while fitting a solar- 
wind like density profile yielded a density of 5xl 2 ° 
cm -3 (where m is the harmonic number) at an 
altitude of 1 Rs- 

Furthermore, we speculated that the solar-cycle 
averaged radiated Type III energy output could 
be similar to that of radiated non-thermal bremm- 
strahlung from accelerated nanoflare electrons. 
We further suggest that escaping electron beams 
could be a viable mechanism for carrying energy 
away from the corona, in a quantity similar to 
the energy introduced into it through nanoflares, 
though this result depends heavily on as yet ill- 
known radio radiative efficiencies. 

It is likely that statistical testing of the differ- 
ent theoretical models of the gen eration and prop- 
agation of Ty pe I II bursts (e.g. bulk et al.l jl979t 



iy pe J 

Duncanlfl979t iRobinson fc Cairnslll998l ) could 



be done. In particular refractive effects on ap- 
parent source sizes, locations, and fluxes. Possi- 
ble other future work includes the usage of sub- 
second or the full 128 ms resolution data, if it ever 
becomes easily usable (e.g. stored on a file sys- 
tem instead of on tapes). After the development 
of a good (and automated) Type III burst detec- 
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tion and discrimination (e.g. Lobzin et al. 2009), 
this work could be repeated almost verbatim. In 
the meantime, it is envisaged to do a statisti- 
cal comparison between the Type III bursts in 
this work and the X-ray flare parameters found 
in the RHESSI flare or microflare lists, and get 
observational constraints on radiative efficiencies 
by comparing with in-situ measurements of elec- 
tron beams associated with the metric/decimetric 
Type III emissions presented in this work. 
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A. Convolution of two elliptical gaussians with different tilt angles. 

The convolution of a gaussian elliptical source (semi- major axis a trU e, semi- minor axis b trU e, and tilt angle 
@true) with a gaussian elliptical beam (abeam, bbearm dbeam) yields an observed source that is also an elliptical 
gaussian (a b s , b b s , obs ). Through Fourier transformation, the following can be derived: 

( trUe) («L - b lbs) cos(2 6 obs ) (al am - bl a J cos(2 6 heam ) {Al) 

a true = \ l(a 2 o bs + b2 obs- a lea m - b lea m + A ) ( A2 ) 



b tme = \l\{a 2 obs + b 2 obs -al eam -bl eam - A) (A3) 



with: 



A _ Kbs - b lbs) cos(2fl obs ) - (a 2 beam - b\ eam ) cos(2g fce 

COs(2 Otrue) 



b obs 



(A4) 



or A _ ( Q L - b lbs) sin(2g obs ) - (a 2 beam - b 2 beam ) sin(2 9 beam ) 

sm(2 6 true ) 

(use whichever has non-zero denominator). 

Moreover, an important result, which can be derived by combining Eqs (A2) and (A3): 



— 3 true + s beam (A6) 



where sf = a 2 + b 2 , with i being obs, true, or beam. I.e., when dealing with the convolution of elliptical 
gaussians, the rms size of the observed source is simply the quadratic mean of the rms size of the beam and 
of the true source, independent of any tilt angle between the true source and the beam. 



B. Non-thermal bremsstrahlung efficiency in the case of the thick-target model: 

For an injected electron spectrum with negative spectral index <5 and low-energy cutoff E Cl and assuming 
non-relativistic cross-sections for bremsstrahlung (Bcthc-Heitlcr) and energy losses, the ratio of total emitted 
X-ray energy over accelerated electron energy can be approximated by the following analytical expression, 
derived by considering the photon spectrum to be a power-law above a low-energy cutoff E Cl and a flat 
spectrum below it: 

4 q ? 1 B{6 /2,1ft) M + 



""v,r'^^j^ f + { 1+ t^) (B1 » 

where a is the fine structure constant (« 1/137), A is the Coulomb logarithm («20-30 for fully ionized 
coronal plasma), z 2 is the average coronal atomic number-squared (rsI.44), m e c 2 =511 keV, B is the Beta 
function, 8 is the electron flux power-law negative spectral index, E c is the low-energy cutoff expressed in 
keV, and / is a slowly-varying function that depends mostly on 5, and is equal to about 0.3 for <5=4. For 
example, for <5=4 and E c —10 keV, one gets a bremsstrahlung emission efficiency of 8.5 x 10~ 6 . We have found 
good (within a 10-20%) correspondence with (exact) numerical computations. 



C. Magnetic field tilt estimation through simple geometry arguments 

From the observation that radio sources seem to be displaced westward by ~2' on the average, we derive 
an average tilt angle of the magnetic field (or at least the direction of the main Type III emission) with the 
local radial to Sun center. 



1G 
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EARTH 



Fig. 17. — Type III emission ge- 
ometry. 



For simplicity, we consider in the following the geometry to be 
wholly within the ecliptic plane, and the Type III emission to have 
infinite directivity (zero beamwidth). 

The "average" source is observed to be at a westward angle \ 
from Sun center (\ ~ 2'), at an altitude Rs + h, where Rs is 
the solar radius. We call the angle between the Sun-Earth line 
and the radial to the source, and a the angle between the radial 
to the Sun and the direction in which the Type III is emitted 
(Figure [TT|) . Then, the following must be true for the emission 
from the "average" source (which is near the solar surface) to be 
observed at Earth: 



<j) + X 



A bit of trigonometry yields: 



sin(co) 



D-(R S + h) cos(0) 
R s + h 



tan(x) 



(CI) 



(C2) 



where D is the Sun-Earth distance («215 Rs)- This is a tran- 
scendental equation which can be iteratively solved for <p. a is 
then simply derived from Eq IC1I The following explicit formula, 
derived using simplified geometry, 



a ~ x + arcsin 



( X-D 



(C3) 



has accuracy better than 1% for h < 0.6Rs- 
Using x=2', and for h between 0.1 and 0.5 Rs (approximate altitude up to which decimetric emission can 

be found, cf. Figured]), one obtains a between 6.6 and 4.8 degrees. Say a=6±l degrees. 

These results were derived from a very simple 2D model, assuming the Type III emission had infinite 

directivity, but they are in good agreement with the Monte-Carlo approach presented in Appendix IDl 



D. Magnetic field tilt estimation through a simple Monte-Carlo approach 

In addition to the simple geometrical approach described in Appendix [Cl we ran a simple 3D simulation to 
estimate the effects of a tilt angle a of the radio emission with respect to the local radial, the half- width Act 
of the cone of emission, and the source altitude h. A hundred thousand sources were uniformally distributed 
in longitude, but only within 15 and 30 degrees in latitude (which reflects reality). For each source, if Earth 
(and the observer) was within the cone of emission, the apparent position of the source on the Sun was 
recorded. If not, it was omitted. In Figure [151 we display the variations in averaged source East- West 
position, if one varies a, Aa, or h. 

Using canonical values of a— 6 degrees, Aa=30 degrees, and h=0.5 Rs, one obtains a westard displacement 
of 137", i.e. very close to our reported observations. It can be said from Figure fTHl that the result does not 
strongly depend on the exact value of h, i.e. a change of 0.1 Rs leads to a change of the average displacement 
by «10". Reasonable (<80 degrees) values of Aa lead to variations of up to ±25" at the most. 

We conclude that the 6±1 degrees tilt angle obtained in Appendix [Cl is further validated by the simple 
3-D model presented in this appendix. 
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Tilt angle [degrees] Half-cone angle [degrees] Source height [R SUN ] 

Fig. 18. — Westward displacement of the average source position, obtained through a Monte-Carlo simula- 
tion. See text for details. 

E. Rank-order correlations 

Table [5] displays the Spearman Rank-order cross-correlation coefficients between burst parameters, for 
bursts observed at 164, 327, and 432 MHz. Spearman's Rank-order cross-correlation acts on the ranking of 
the values instead of on their actual values, as a regular cross-correlation (e.g. Pearson's) would. This has 
the huge advantage that even a non- linear (polynomial, exponential, power-law, etc.) relationship between 
parameter x and parameter y can achieve a cross-correlation coefficient near unity. The parameters are: 

• r is the angular distance to Sun center. 

• Ts.true is the peak brightness temperature (averaged over 10 s) of the deconvolved gaussian source (as 
described in Figure ITUl and Section l3~T]) . 

• S v is the burst flux (averaged over 10 s). 

• Stme is the FWHM size of the source, deconvolved as explained in Section |3~T1 

• Otme is the orientation of the semi-major axis of the source, from the map x-axis. 

• (a/b)true is the ratio of the semi-major to semi-minor axis of the source, an easier quantity to visualize 
than the eccentricity. 

/stme is a measure of the source's radial elongation, normalized to source rms size. 

• Strue.azirnuthal 

I stme is a measure of the source's azimuthal elongation, normalized to source rms size. 
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